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We study the peculiar velocities of density peaks in the presence of primordial non-Gaussianity. 
Rare, high density peaks in the initial density field can be identified with tracers such as galaxies 
and clusters in the evolved matter distribution. The distribution of relative velocities of peaks 
is derived in the large-scale limit using two different approaches based on a local biasing scheme. 
Both approaches agree, and show that halos still stream with the dark matter locally as well as 
statistically, i.e. they do not acquire a velocity bias. Nonetheless, even a moderate degree of (not 
necessarily local) non-Gaussianity induces a significant skewness (~ 0.1 — 0.2) in the relative velocity 
distribution, making it a potentially interesting probe of non-Gaussianity on intermediate to large 
scales. We also study two-point correlations in redshift-space. The well-known Kaiser formula is still 
a good approximation on large scales, if the Gaussian halo bias is replaced with its (scale-dependent) 
non-Gaussian generalization. However, there are additional terms not encompassed by this simple 
formula which become relevant on smaller scales (k > 0.01 h/Mpc). Depending on the allowed level 
of non-Gaussianity, these could be of relevance for future large spectroscopic surveys. 
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I. INTRODUCTION 

Observations of the large scale structure in the Uni- 
verse that use galaxies, clusters or other tracers of the 
density field are done in redshift space: the distance 
is generahy inferred using the redshift z, which receives 
a contribution from the line-of-sight velocity of the ob- 
ject. These velocities arc due to the gravitational field 
which is correlated with the density field itself. On large 
scales where linear perturbation theory in the density 
field applies, the leading contribution is the squashing 
(or stretching, in case of an underdensity) of a volume 
element in redshift space relative to real space. In this 
limit, there is a simple relation between the real- and 
redshift-space power spectra, Pg and Pg,s, respectively, 
of a tracer 'g' [i|, Q: 



Pg,.(fc,M)=(l + £M-^) Pg{k), 



(1) 



where ^ is the cosine of the k vector with the line of sight, 
/ = d In D/d In a is the logarithmic derivative of the lin- 
ear growth factor, and 6i is the linear bias of the tracer 
population. Apart from the large-scale, small-correlation 
limit, the relation Eq. ([1]) makes two assumptions: first, 
that the tracer population is characterized by a deter- 
ministic local bias. In particular, if the tracer density 
(5g(x) is a local function F{6{x.)) of the matter density 
perturbation 6, we can expand in (5(x) to obtain [3| 



'2 £-2/ 



Sg{x.)^biS{x) + ^6\x.) + . 



(2) 



where the bias parameters are either to be seen as free 
empirical parameters, or can be determined using vari- 
ous theoretical approaches. Local biasing amounts to the 
assumption that dark matter halos form in high-density 



regions (peaks) in the initial density field. This assump- 
tion holds well in the high-peak / massive halo regime, 
which we assume throughout. Hence, in the following 
we will somewhat loosely use "peaks" and "halos" inter- 
changeably. 

The second assumption used for Eq. ([T]) is that the 
cosmological density field is Gaussian on large scales. 
While we will retain the local biasing scheme, we are 
interested in relaxing the second assumption of Gaussian- 
ity. Recently, there has been renewed interest in probing 
the Gaussianity of the primordial seed perturbations via 
large scale structure (see Q for a review). The simplest 
way of obtaining a non-Gaussian field is by adding a local 
non-linearity: 



$(x) = $G(x) + /NL(<i>^(x)-($^)), 



(3) 



where $g is a Gaussian random field, and $ is the re- 
sulting non-Gaussian field. Of course, one can add higher 
powers to the scries Eq. ([3]), though the quadratic term 
usually has the largest impact. Following standard con- 
vention, we let $ stand for the primordial potential, re- 
lated to the density field through the transfer function 
and Poisson equation (see App. E)) . 

As shown by Q and confirmed by [a, 0] i in the presence 
of non-Gaussian initial conditions of the local type, halos 
acquire a scale-dependent correction to their bias which 
becomes important on large scales: 



6i^fei+2/NL(6i-l)<5cA^(fc), 



(4) 



where M{k) oc fc~^ is the relation between density and 
potential in Fourier space (see App. |A| . 

Besides the local model, several other possible bispec- 
trum configurations have beenproposed in the literature, 
such as the equilateral (e.g., Q) and folded types (e.g., 
Q). We summarize these types of non-Gaussianity and 



the relation between the potential and matter perturba- 
tions in App. 1^ 

Given the significant impact of (local) non-Gaussianity 
on the power spectrum of biased tracers, it is then natural 
to ask what happens to the power spectrum in redshift 
space, and whether Eq. ([T]) still holds. Furthermore, the 
distribution of relative velocities between tracers at Xi 
and X2, 



Su = u(x2) - u(xi) 



(5) 



is itself an interesting probe of non-Gaussianity [10|. 

In this paper, we assume sub-horizon scales through- 
out, and adopt the Newtonian gauge; further, we work 
in Lagrangian coordinates. At first order, the transition 
to Eulerian coordinates (which all observations as well as 
simulation measurements are made in) simply amounts 
to replacing the Lagrangian bias parameters with their 
Eulerian counterparts, in case of the linear bias simply 
bE,i = bL.i + 1- 

While we will focus on the local type of primordial 
non-Gaussianity for the most part, the expressions ob- 
tained can easily be evaluated for any given primordial 
bispectrum, and we will show selected results for other 
bispectrum shapes. 

The paper is structured as follows: ?[II1 introduces den- 
sity and velocity correlations, and some notation. § IIIII 
contains two different derivations of the moments of the 
relative peak velocity distribution. We discuss the distri- 
bution of peak velocities in § IIVI Finally, § |V] presents 
the power spectrum of peaks (halos) in redshift space. 
We conclude in 5 IVll 



II. PRELIMINARIES: DENSITY AND 
VELOCITY FIELDS 

We consider the linear ovcrdcnsity and velocity fields 
of matter <5(x), v(x) as random fields, related by the (lin- 
ear) continuity equation. 



(5 + Vv = 0. 



(6) 



Since for the linear density field S = aHf S, we can define 
a scaled velocity u which satisfies 



u ; 



aHf 



^ Vu = -S. 



(7) 



In the following, we will always deal with smoothed den- 
sity and velocity fields which are indicated by a subscript 
R: 

Sni^) = |dVw^i?(|x-y|)<5(x) (8) 

Uii(x) = yd-VW^i?(|x-y|)u(x), (9) 

where Wn is a normalized window function. Specifically, 
we use a real-space tophat filter for Wr, though the shape 



of the window function has negligible impact on our re- 
sults. Following the standard convention, we choose the 
smoothing scale R corresponding to a halo of mass M 
to be determined hy R = (SM/Anp)^/^, where p is the 
background matter density. 

Below, we will need two-point correlations of the den- 
sity, velocity, and the cross-correlation between the two. 
Using the smoothed matter power spectrum. 



Pfl(fc) = W^k) P{k) 



(10) 



where P{k) is the unsmoothed power spectrum and W is 
the Fourier transform of the window function, the two- 
point correlations are given by: 

^fl=^f(r) ^ (5fi(xi)<5fl(x2)) (11) 

= ^ldkk'Pn{k)jo{kr) 

C^"(r) ^ {5n{^,)f-nn{^2)) (12) 



1 

2^' 



2 , dkkPn{k)ji{kr) 



CrM ^ -(Ufl(xi).Ufl,(x2)) 



(13) 



-^ j AkPn{k)3o{kr). 



Here, r = |x2 — xi| and f = (x2 — xi)/r. S,j^{r) denotes 
the cross-correlation between the density at point xi and 
the matter velocity at point r2 projected along the sep- 
aration vector. It is negative, since in the presence of an 
overdensity (5fl(xi) > the streaming motion will be di- 
rected towards point xi. Finally, we define the variance 
of the smoothed density field a^ and the one-dimensional 
smoothed velocity dispersion cr„: 



4^^«(0); 



2 t^-uu 



a"(o). 



(14) 



We will also need expressions for various three-point cor- 
relators of the line-of-sight velocity and density. The 
density three-point function is given in terms of the 
smoothed matter bispectrum Bj^, which is defined in 
App. El by 



cnxr,x2,x3)=y^y^ 



d^k^ 

r)3 



(15) 



„ifk 



i.(xi-X3)+k..(x.-X3)] 5^(ki,k2,-ki -k2). 



Since we are interested in two-point peak correlations, 
we will always encounter degenerate triangles with X3 = 
X2. If instead two other vertices coincide, we can always 
relabel the indices to bring the correlation into this form. 
We define the general, mixed density/ velocity degenerate 
three-point correlation as 






[r) = 

X«(x) = 



(aW(xi)xM(xi)X("Hx2)), 
^fl(x), 1 = 

f •Uj^(x), I = 1, 



(16) 
(17) 



where r = X2 — Xi. Using that in Fourier space, ur ■ f = 
k-r / {—ik)5fi, and noting our definition of r, we can write 
these correlations as 



crir) 



d-'k^.kr [ fe-^ 



(27r)3 



d''k2 I k2 ■ r 



(2^) 



-ife 



-(k + k2)-f 



-ilk- 



, ,2 ; Bfl,(k + k2,k2,k). (18) 

k2 r 



The corresponding Fourier-space quantity can be defined 
by 



er{r) = 



Imn 



(fc) = 



d^k 

k) J (2^)3 

k2z 

k. 



-^ik-r ■l-\-m-\-n hnn 



(k) 



|k + k2P 



(19) 



^) B;7(k + k2,k2,k), 



where here and throughout we use the small-angle (flat 
sky) approximation, and let the line-of-sight direction be 
along the z-axis. Further, we have taken out powers of i 
in order to make j^'™" real. 

For the numerical results, we adopt a flat ACDM cos- 
mology, with 17™ = 1 - r^A =: 0.28, VLb ^ 0.046, crs = 0.8, 
and Hs = 0.958. The Gaussian Lagrangian bias of ha- 
los is calculated as h^ = 5c/fy\, where 5c = 1.686 is the 
linear collapse threshold. Note that our expressions are 
more general in the sense that they can be used with any 
local bias parameters (see § IIIIAp . 



III. VELOCITIES OF DENSITY PEAKS 



velocity with respect to the large-scale flow of matter. 
In Gaussian N-body simulations, this ansatz has been 
shown to hold well on large scales, and one would not 
expect that mildly non-Gaussian initial conditions will 
change this. 

On the other hand, as shown in [illll2|, density peaks 
have a statistical^ scale-dependent velocity bias when the 
peak constraint is employed. Here, statistical means that 
each halo still moves with the dark matter locally (no 
local velocity bias), but a velocity bias still comes about 
from the fact that peaks in the density field arc in special 
locations. This velocity bias from the peak constraint 
scales as fc^, and we will neglect it here for simplicity 
since we are interested in large scales. However, it is not 
apparent whether a non-Gaussianity of the density field 
on large scales can also source such a statistical velocity 
bias. We defer this question to § IIIIBl and assume in 
this section that there is no such bias. 

Using the ansatz of unbiased velocities, it is easily pos- 
sible to derive the first few moments of the velocity differ- 
ence distribution. For this, we expand the peak density 
to second order (in Lagrangian space): 

<Spk(x) = 6i5fl(x) + ifeL,2<5|(x), (21) 



bl 



where the second order bias parameter is hL,2 = OJ, m 
the high-peak (high-significance) limit we work in (see 
below for the general result). For the sake of clarity, we 
will omit the subscript R for the smoothed fields, and let 
5i ~ 5fi{xi) and similarly for u for the remainder of this 
section. 

By assumption, Upi^ = u, and we have for the relative 
velocity (to first order in the correlations): 



In this section we derive the moments of the distribu- 
tion of velocities of density peaks in two different ways, 
both based on the local bias model of Eq. ^ (i.e., local 
in terms of the non- Gaussian density field S). We are 
interested in the relative velocity of two peaks projected 
on the separation vector. 



[U(x2)-U(xi)] -f. 



(20) 



This is the relevant quantity to compare with observa- 
tions of large scale velocities, in particular redshift-space 
distortions (see § |V| . Homogeneity and isotropy dictate 
that the statistical properties of 6u such as its moments 
only depend on r = |x2 — xi|. 



A. Unbiased velocity ansatz 

The first, much simpler approach is to assume that the 
velocities of peaks are locally and statistically unbiased. 
The absence of a local velocity bias is motivated by the 
fact that local processes such as gravitational collapse of 
small-scale ovcrdcnsities should not induce any relative 



{Su)ir) 




bLSl + -bL,25!][l 



+ hLh + 7^hL,25l] 



X f • [U2 - Ui] 



26ier(xi,X2,X2) 



26Leg'(0 

+ &L,2Cfl"(xl,Xi,X2) 

2bU'R{r) - 26ie^"*>) + hierir)- (22) 



I 2 t-SSui 



Here we have used Eq. (|16p to express the various three- 
point correlations in a uniform way. Note that when 
reordering the arguments of the three-point function in 
order to bring them into the form of Eq. ()16p . this can 
entail a change of sign when swapping xi o X2 for terms 
odd in u (due to the projection onto f). 

We see that while we recover the standard expression 
in the Gaussian case (obtained by setting all three-point 
functions to zero; see e.g. [li|), there are several three- 
point terms contributing to the mean relative velocity. 
Note however that {5u) is non-zero only for biased tracers 
^L 7^ even in the presence of non-Gaussianity. 



Next, the dispersion of velocity differences is given by 

{6u^){r) = ([l + 6i,5i][l + 6i(52](f.[u2-ui])2) 

= 2<jI + 2er(r) + 26i(5i(f • [u2 - Ui])^) 

+2bL[CR\0) + er'(0 - 2^fl""(0]- (23) 

Again, we recover tlie standard expression in tlie Gaus- 
sian case, with several additional three-point contribu- 
tions. These contributions vanish however for an unbi- 
ased tracer (6l = 0). 

It is straightforward to generalize {6u) and (Su^) to 
the case when two different tracers with arbitrary linear 
and quadratic bias parameters are considered. The first 
two moments of the relative velocity distribution between 
different tracers are then given by: 



{Su)ir) = {ha + bib) en (r) - 26i.6ibCrW (24) 

{Su'){r)=2al + 2Cn''ir) (25) 

+ iha + bit) Kr'(O) + Cn'ir) - 2ei""(r)]. 

Here, 6ia,&ib denote the linear (Lagrangian or Eulerian) 
bias parameters for the two tracers, while 62a, &2b denote 
the corresponding quadratic biases [Eq. ^]. 

Finally, there is also a third moment of Su, resulting in 
a skewed velocity difference distribution as already no- 
ticed by [lOl ■ It is given by 

{Su^){r) = ([l + 6L(5i][l + 6i<52](f-[u2-ui])3) 

= ([f.(u2-ui)]3)=6er"(0, (26) 

where for the second equality we have again assumed the 
absence of any four- and higher-point terms. The third 
moment of velocities is obviously absent in the Gaussian 
case, but for non-Gaussian initial conditions it is non-zero 
even for unbiased tracers (e.g., matter itself). In the next 
section we outline the second derivation of these results 
(more details can be found in App. IB)) . 

B. Derivation in statistical field theory 

The assumption that peak velocities are still sta- 
tistically unbiased in the presence of large-scale non- 
Gaussianity seems natural. On the other hand, identi- 
fying peaks in the density field is a non-linear process, 
and thus it is desirable to have a proof for this assump- 
tion. Our second derivation of the distribution of rel- 
ative velocities provides such a proof without making 
any assumptions apart from the large-scale limit (i.e., 
small correlations), and local biasing in the physical, non- 
Gaussian density field. 

In accordance with local biasing (see e.g. [3, [i3])j we 
associate peaks with positions where the smoothed den- 
sity field dji is larger than a threshold vcfr (hence, for 



1/ = 3 we are selecting 3a peaks of the smoothed density 
field). The average number density of such peaks is then 
given by 



"pk = (6[(5_r(x) - lyo-Rl), 



(27) 



where O is the Heavyside function, and (•) here denotes 
the ensemble average. Similarly, the correlation function 
of peaks is given by (using the same shorthand notation 
as above) 



c /\ {0[Si - iyaR]Q[S2 - v(tb]) ^ 



(28) 



'pk 



We are interested in the pairwise velocity of peaks along 
the separation r. We can write the pair-weighted proba- 
bility distribution for the two velocities Uri, Ur2 as: 

7^(Url,Ur2;Xi,X2) = Tlpj^^ X (29) 

(^^(f • Ui - Url)5D{r ■ U2 - Ur2)Q[Sl - i^Crfl]e[(52 " I^Tb])- 

Note that by construction, the peak correlation function 
Eq. pS)) is recovered through 



1+Cpk(?') = / dUrl / dUr2'P{Ur2,Ur2;Xl,X2)- (30) 

Given Eq. (j^ . the volume- weighted distribution Vsu of 
the pairwise velocity 6u = Ur2 — Uri at a given separation 
is obtained as 

7'5u(u;Xi,X2) = -— — / dUrl'P{Url,U+Url]y.l,y.2)- 

1 + Cpk(r) J 

(31) 
We now write the ensemble average of any functional 
J^[u] of the velocity field as a functional integral over all 
possible realizations of the (linear) velocity field u(x): 



(^[u])=yp[u(x)]7'[u]J-[u] 



(32) 



We make use of a result of statistical field theory (see 
e.g. |13| - K15| for related applications), expressing the par- 
tition function Z[i] in terms of higher order (connected) 
correlation functions: 

Z[J] = /r'[u(x)]7'[u]exp|i /d^x J(x)-u(x)l (33) 

= exp|f;^y'd3xi...y'd3x„ 

X 4".'..„(xi,...x„)J'Mxi).../"(x„)l. (34) 



Here, 



.(n) 



4i...i„(xi'---^«) = (Wji(xi)...Wj„(x„))con (35) 

is the connected n-point correlation function of the un- 
smoothed velocity field. We can now use Sr — ~^vlr 



together with Eq. ((33|) to obtain the velocity probabihty 
distribution of peaks defined in Eq. (P^ . in close anal- 
ogy to the derivation for density correlations presented 
in [13 • This calculation is detailed in App. IB] 

We only keep terms linear in the correlations, an ap- 
proximation valid on large scales. Once we have the ex- 
pression for Psu [Eq. (|5T|) ]. it is straightforward to obtain 
the moments of the relative peak velocity Su via 



10 



(5M")(r) = / dSuSu'^VsuiSu-^r). 



(36) 



Note that the central and non-central moments agree to 
first order in the correlations (since (du) is already first 
order). We then obtain (App. IBJ): 



{Su'' 



{Su)-- 


= '^bLCnir) 

+ &i[ei'"(r)-2ei"^(r)] 


(37) 


{Su') = 


==2[a2+eF(r) + 6Ler'(0) 


(38) 


{Su') = 


-eerir) 


(39) 


SuY = 


= 0. 


(40) 



These moments agree with those derived in § IIII Al Note 
that there is no connected fourth moment (kurtosis) at 
this order, since we have assumed no primordial trispec- 
trum. The fact that the moments of the peak velocity 
distribution agree with those derived in § IIII Al shows 
that no large scale statistical velocity bias of peaks is 
induced by non-Gaussianity. 



IV. DISTRIBUTION OF RELATIVE 
VELOCITIES 

Fig. [1] shows the moments of Su in the Gaussian and 
non-Gaussian case as function of the separation r. Here, 
we have assumed local non-Gaussianity with /nl = 100, 
and halos with 10^' Mq/H at z = 0, corresponding to a 
smoothing scale of i? = 3.1 Mpc/h and a Lagrangian bias 
of bL = 0.88. 

The mean streaming velocity of halos (Su) is always 
negative, since overdense regions have a net infall ve- 
locity. It is clearly enhanced on large scales by non- 
Gaussianity. This is mainly due to the term &|Cfl*", 
which corresponds to the density-weighted analog of ^^" 
and is negative for positive /nl, i-e. for a positively 
skewed density field. The effect on the variance {Su'^) on 
the other hand is very small. In Fig. [TJ we have split the 
variance into a scale-dependent and scale-independent 
part, and denoted the latter by {Su'^)o. The scale- 
independent standard deviation of Su is 7.9 Mpc//i for 
the adopted smoothing scale and /nl = 100, while it 
is 8.1Mpc//i in the Gaussian case, a difference of only 
3%. Thus, the effect of non-Gaussianity on the scale- 
independent as well as scale-dependent part (Fig. [T]) of 
the variance is small. Note that the non-Gaussian terms 
in the variance are proportional to the bias 6^. 
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FIG. 1: Moments of the relative velocity distribution Vsu 
as function of the separation r. The thick lines show the 
moments for a cosmology with primordial non-Gaussianity 
of the local type (/nl = 100), while the thin dotted lines 
show those for vanishing non-Gaussianity. We have separated 
the variance into a scale-independent part (5u^)o and a scale- 
dependent part (see text). In all cases we assumed halos with 
M = 10" Mq/Ii at z = (smoothing scale R f» 3.1Mpc//i, 
bL = 0.88). 



Finally, the third moment of Su is significant. The 
skewness, defined as {Su') / {Su'^)'^'^ , is between 0.1 and 
0.2 on a wide range of scales. The skewness which quan- 
tifies how strongly non-Gaussian the velocity distribution 
is, depends only weakly on the halo mass, through the 
smoothing scale R, and redshift. Since {Su'^) slowly de- 
clines with increasing smoothing scale while {Su') is es- 
sentially independent of R, the skewness is slightly larger 
for higher mass halos, e.g. 10% higher for 10^'' Mq/H ha- 
los compared to 10^'^ Mq/H. Note that for typical choices 
of the smoothing scale, there will be a significant non- 
linear correction to the variance {Su^). 

Further, it is important to note that the skewness is 
non-zero for other bispectrum shapes as well. For exam- 
ple, a bispectrum of the equilateral type (see App. \K\ . 
again with /nl = 100, yields a skewness of about 0.06 
at r = 100Mpc//i, only a factor of three smaller than 
the skewness obtained in the local model. On the other 
hand, the power spectrum of halos only receives a scale- 
independent correction in the equilateral model (jla). 
§|V|). Thus, the distribution of relative velocities could 
be an interesting avenue to test especially those models 
of non-Gaussianity which can only be weakly constrained 
from the clustering of halos. 

Generally, the skewness is positive for positive ^™" (?-) ^ 
which is usually the case when S,^^^ is positive (i.e., for 
positive /nl)- At first this might seem counterintuitive. 
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FIG. 2: Distribution of the relative velocity Su for a separa- 
tion r — 100 Mpc/h in the non-Gaussian (solid) and Gaussian 
(dotted) case, for the same parameters and halos as in Fig. [l] 
Note that an equilateral model with /nl ~ 300 would lead to 
a result very similar to the local model shown. 



since for a positively skewed density field there are more 
overdense than underdense regions. Hence, one expects 
that velocities projected on f are preferentially "inward" , 
i.e. negative. However, this is only true for density- 
weighted velocity correlations, such as ^^^" which is in- 
deed negative. ^™" = {u\^U2r) on the other hand is 



a velocity-weighted velocity correlation. ^™" > thus 
means that high velocity regions on average move apart 
(while there is no such net motion in the Gaussian case). 
High- velocity regions correspond to infall regions of over- 
densities, so ^^"^ corresponds to enhanced infall motions. 
This is not unexpected for a density field with more high- 
density regions. Of course, the converse holds for a neg- 
atively skewed density field (/nl < 0). 

It is also interesting to look at the full distribution of 
'Psui^u). Since by assumption the fourth and higher ve- 
locity moments are small (and in fact all moments other 
than the variance are small on large scales), we can em- 
ploy the Edgeworth expansion (e.g., [1| 



l es) , w 
,113): 



P5«('5'«;xi,x2) = <^ 1 



{Su)Su {Su^ 



2cjI 



8af, 
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{6z 



12cr3 4au 
1 






-6 



V47ro^ 



:exp 



3W 

4^2 



7T7 • (41) 



We also arrive at this expression directly using the sta- 
tistical approach of App. [B] The distribution of the rel- 
ative velocity Su at a separation of 100 Mpc//i is shown 



in Fig. [2] in the Gaussian and non-Gaussian case, for the 
same halos and local /nl = 100 as in Fig. [TJ The ef- 
fect of non-Gaussianity is clearly visible as skewing the 
distribution towards positive velocities, as expected from 
the positive third moment. Correspondingly, a negative 
/nl would lead to a velocity distribution skewed towards 
negative relative velocities. Again, these results arc only 
weakly dependent on halo mass and redshift. Further- 
more, the relative velocity distribution of the same halos 
in an equilateral model with /nl ~ 300 looks very similar 
to the local model result shown in Fig. [21 

Recently, observational studies have found evidence for 
larger bulk flow motions than expected in ACDM. These 
studies have used galaxy peculiar velocity surveys p^ . 
the kinetic Sunyaev-Zeldovich (SZ) effect [I3, and the 
estimated initial relative velocity of the merging halos in 
the Bullet cluster [201 ■ l^ig- [D shows that even a mod- 
erate primordial non-Gaussianity that is still allowed at 
the 2(7 level can signiflcantly alter the tails of the veloc- 
ity distribution. Thus an obvious question is whether the 
existence of primordial non-Gaussianity can alleviate the 
tension between the reported bulk flow observations and 
the ACDM scenario. Unfortunately, it is not straightfor- 
ward to relate our predictions of linear perturbation the- 
ory in the large-scale limit to these observational results. 
Thus, we do not attempt any quantitative comparison 
here. 



V. MATTER AND PEAK CORRELATIONS IN 
REDSHIFT SPACE 

One of the most important observational effects of 
large-scale velocities is their impact on correlations of 
tracers such as galaxies, clusters, or the Lyman-a forest, 
through redshift distortions. Positions are measured in 
redshift space, whose coordinates s are related to real- 
space coordinates x by a shift along the line-of-sight di- 
rection: 



S = X + fUgZ. 



(42) 



Here, we have chosen the line-of-sight to be along the 
2;-axis. Throughout, we will work in the flat sky limit. 
The statistics of the density fleld (of matter or peaks) 
in redshift space can be determined using the real-space 
statistics of the density and velocity flelds together with 
the mapping Eq. P^ . 

We start with the normalized distribution Psu{Su,r). 
The correlation function in redshift space ^^ can then be 
written as a convolution by the velocity-difference distri- 
bution 'Psui6u;r) as |2l.l2]| 

1 + Usz,s^) =. fdSu[l + ^{riSu))]rsuiSu;r),iA3) 



r{Su) = ,Js\ + (s. - f^fSu)'^, 



where the factor fi = Vz/r multiplies Su, the total velocity 
along r, to obtain the velocity along the line-of-sight. 



This expression is exact in the flat sky hmit. On large 
scales, the displacements from real to redshift space are 
much smaller than the separation r, and we can expand ^ 
as well as Vsu around r^ = Sz- Generalizing expressions 
given in e.g. (2l|, we obtain: 



1+^s(Sz,S_l) 



E 

n,7n—0 

d 



(-M/) 



71 + 771 



7 ' 7T) ' 



(44) 



X 



X^(1 + ^(^))X^('5""^™)(^)' 



ds 



dsl 



where now /i = s^/ s. Further, we keep only those terms 
which are first order in large-scale correlations. This gives 



6(s.,sJ = i{s)-f^i{5u)'{s) + ^{5uY{s) 



3, ,3 



/v. 

6 
6 



{5v: 



3\iii I 



s) + 



2 



{5u^){s)e{s) 



{Su'){s)C"is), 



(45) 



where primes denote derivatives with respect to s^. The 
last two terms are formally second order, however the 
moments of the velocity {Su"^) can be significant even on 
large scales due to small-scale motion. However, we ver- 
ified that the last, non- Gaussian term is negligible com- 
pared to the other terms. Since the Gaussian quadratic 
term is quite small on large scales as well, we will neglect 
those two terms in the following. 



We can now Fourier transform Eq. (j45j) term by term, 
noting that every derivative with respect to Sz brings 
down a power of ikz = ifik. Defining the Fourier trans- 
form of the nth velocity moment (in the flat sky limit) 



{Sin{r) = 



d^k 

(2^ 



((5u")(fc), 



(46) 



we have: 



P{k) - fin^k{6u){k) + i-{-fj,'ik^){5u^k)) 



6 



■^i^{ikf{Su^){k). 



(47) 



This expression applies to both peaks and matter. With 
the results of § IIIIl it is then straightforward to write 
down the redshift-space peak power spectrum Eq. (|T7)) 
in the presence of non-Gaussianities. Again, these re- 
sults apply to bispectra of arbitrary shape- and scale- 
dependence. 

Fourier transforming the 6u moments given by 
Eq. (USD, Eq. ^, and Eq. ^ is easily done by re- 
placing the correlations ^'™"(r) with ;>''"" (fc). The full 
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FIG. 3: The terms of Eq. (|52|) which contribute to the 
redshift-space power spectrum and are not captured by the 
Kaiser formula, as function of k. We also show the lead- 
ing real-space term p (k) for comparison. All terms are for 
10^^Mq//i halos at z = and local non-Gaussianity with 
/nl = 100. We have set /i = 1 and divided all terms by the 
matter power spectrum. Note that non-linear corrections to 
these linear predictions are expected above k '--^ 0.05 h/Mpc. 



expression of the large-scale limit then reads 



~fti^ik[2bLiik)-^PRik) 

+foi(z/*"(fc)-2»/"^(fc))] 

-^/x'A;2[2(ifc)-2p^(fc) 



26l;>""'(A) 



6 



-f4foip^""(fc)] 
■fi^{ik)^i^p''''''{k). 



(48) 



The first line contains the real-space peak power spec- 
trum including the well-known non-Gaussian contribu- 
tion (111, see Eq. (|B30)) in App.|B]). Note that in Eulerian 



space, b^ is to be replaced by 6162. Using p 
which follows from Eq. (|19p . we get 



hnu 



IrnS 



Ik, 



Ps{k,fi) = blPnik) + 2bLfti^PR{k) + ffi^Pnik) 
+6i/^^(fc) + V6i/^^(fc) 
~2f„Hlkp'^\k) 
+fti%Lk\''^\k)-2p'^'\k)] 
+/V^ k'^p'''"'{k) (49) 

The terms in the first line constitute the usual Gaus- 
sian power spectrum of biased peaks in redshift space. 
Following convention [5|-|7[ , we can write the first term in 
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(see App.lA]). 



the second line as a non-Gaussian correction A5l to the 
Gaussian halo bias (see e.g. Q): 



A55 



(fc) 



6l 2 ^ P,i(A:) 



(50) 



Thus, if we include this no n- Gaussian bias and neglect 
the terms in the last three lines above, the Kaiser formula 
Eq. (m is still valid in the presence of large-scale non- 
Gaussianity: 



P. 



Kaiser 



(fc,^) 



/A^' 



bL + AbL 



(bL+AbLfPRik). 



(51) 

However, there are additional terms which we will denote 
as Pf G. 

pf«(fc,/i) = -2k/-\k)[f^,'bl + .f^,HL] 

+fcV""'W[/V&L + ./V]- (52) 

Corrections of order /nl to the Kaiser formula were 
also found by [221 for the variance of the density field 
in redshift space. Fig. |3] shows the remaining terms 
which violate the Kaiser formula in the presence of non- 
Gaussianity on large scales. Comparing them to the real- 
space contribution p^^^ , it is clear that these terms are 
far subdominant. However, they become increasingly im- 
portant towards smaller scales. In order to fully evalu- 
ate their importance on scales k > 0.05 /i/Mpc, it will 
be necessary to extend the perturbative calculation to 
higher order or to measure the effect in N-body simula- 
tions. In the adopted cosmology, the non-linear scale for 
matter is fcni w 0.27/i/Mpc, where kn\ is defined through 
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FIG. 5: Ratio of the redshift-space power spectrum with pri- 
mordial non-Gaussianity (/nl = 100) to that with Gaussian 
initial conditions as function oi fi at k — 0.1/i/Mpc. Again, 
results are shown for lO'^^ Mq/H halos at z — 0. Shown is the 
actual ratio including all terms (Eq. (|51|l and Eq. (|52p '). and 
the ratio one obtains from the Kaiser formula (Eq. (|51|) only). 
Note that this is mainly illustrative, since we are showing lin- 
ear predictions at a scale already in the quasi-linear regime 
(non- linear correction from halof it ~ 5%). 



A2(fcni) = fcj^iP(fcni)/27r2 ^ ^ ]v|q^p ^j^^^^ ^^^^ ^jjj ^^ somc- 
what smaller for massive halos whose power spectrum is 
enhanced by b^. 

Fig. m shows the same results for the equilateral and 
folded bispectrum shapes. The overall corrections to the 
Gaussian power spectrum are much smaller, since these 
bispectrum shapes are suppressed in the squeezed limit. 
The relative importance of the additional terms from 
Eq. (|52|) is qualitatively similar to the local case. Note 
that for all bispectra, the relative significance of these 
terms compared with the real-space contribution p^^^ is 
independent of /nl at leading order. 

As an example illustrating the effect on the redshift- 
space power spectrum. Fig. [5] shows Ps(fc,/i) at fc = 
0.1/i/Mpc as function of fi, when using the Kaiser for- 
mula Eq. ([STjl and when including all terms. At these 
small scales, the difference is clearly noticeable. This is 
only a rough illustration however, since in reality non- 
linear corrections will contribute significantly at that 
scale: the non-linear power spectrum (from halof it |23| ) 
fc = 0.1 h/Mpc differs by -5% from the hnear P(fc). Stih, 
we expect that the corrections to the Kaiser formula will 
remain at least at the percent- level for k > 0.1/i/Mpc. 
That would make them possibly relevant for upcoming 
spectroscopic surveys such as BOSS and HETDEX, and 
certainly for future space missions such as JDEM and 
EUCLID, which are expected to constrain redshift dis- 



tortions at the sub-percent level [2J, |25| . 

Again, these results do not depend sensitively on halo 
mass and redshift: if, for example, we increase the halo 
mass, there is a partial cancelation between the larger 
smoothing scale and the increased halo bias, leaving the 
relative size of the terms in Eq. (|52p roughly unchanged. 
For z > 0, / increases (asymptoting to 1 at very high 
rcdshifts), but this is compensated by larger biases at 
fixed mass. 

Finally, we note that the last term in Eq. ([5^. which 
comes from the skewness of the velocity distribution, also 
contributes to the redshift-space power spectrum of an 
unbiased tracer, for example matter itself (see also [22l|). 
This is in contrast to the linear matter power spectrum in 
real space, which only receives a much smaller four-point 
correction of order /^l- 

The reasons that the additional terms are so much 
smaller are first, that the dominant effect on the mean 
streaming of halos is taken into account by the mixed 
density- velocity part of the Kaiser formula (when includ- 
ing the scale-dependent bias correction); and second, that 
the third moment of Su, though significantly changing 
the distribution of Su, only affects the correlation func- 
tion via three derivatives [Eq. (j45|) ]. or three powers of k 
in case of the power spectrum. This strongly suppresses 
its contribution on large scales. Note that this is spe- 
cific to the power spectrum, and a larger effect of ^^^ is 
expected e.g. for the bispectrum in redshift space. 

VI. CONCLUSIONS 

In this paper, we have studied the large-scale motions 
of peaks in the density field in the presence of primor- 
dial non-Gaussianity. In the high-peak regime we are 
interested in, peaks can be identified with dark matter 
halos and visible structures such as galaxies and clusters. 
We have derived the statistics of the relative velocity of 
halo pairs in two different approaches based on a local 
biasing scheme. Note that the bias is local in the phys- 
ical, non- Gaussian density field. The first, simpler and 
more accessible approach assumes that on large scales, 
halo velocities follow those of the total matter. Given 
linear and quadratic bias parameters, it is then straight- 
forward to calculate the moments of the relative veloc- 
ity distribution. In the second approach (see App. IBJ). 
no assumptions are made apart from local biasing and 
the large-scale, small-correlation limit. Both approaches 
agree, showing that the assumption of statistically unbi- 
ased halo velocities is consistent on large scales even in 
the presence of non-Gaussianity. 

Interestingly, the presence of non-zero three-point cor- 
relations on large scales leads to significant changes in 
the velocity distribution of halos, and matter itself. For 
a positive /nl, the mean streaming velocity is signifi- 
cantly enhanced on large scales, and a non-zero third 
moment leads to a positive skewness of the velocity dif- 
ference distribution. In the local model, the skewness 
(Su^) I {Su'^f/'^ is of order (0.1 - 0.2) x (/nl/100) on a 



wide range of scales. Even in the equilateral model, which 
has little effect on the halo power spectrum, the skewness 
still reaches 0.06 for the same /nl- This suggests that the 
velocity difference distribution can serve as an interest- 
ing probe of non-Gaussianity (as already pointed out by 
[l^). Furthermore, these findings could be of relevance 
to recent observational reports of significantly higher ve- 
locities than expected in the ACDM (or more generally. 
General Relativity -I- Dark Energy) scenario [18l - l20 |. In 
order to evaluate this quantitatively however, one has to 
take into account non-linear corrections, via perturbation 
theory and/or simulations. 

Our second result is an expression for the redshift- 
space power spectrum of halos (or matter) in the presence 
of non-Gaussianity. The well-known Kaiser formula re- 
lating real- and redshift-space power spectra, extended 
by the scale-dependent halo bias (see e.g. Q), receives 
non-Gaussian corrections which become relevant on small 
scales. Note that through these corrections, the redshift- 
space power spectrum measured at a given scale as func- 
tion of linc-of-sight angle in principle allows for a direct 
measurement of non-Gaussianity from a given tracer pop- 
ulation, without any reference to the underlying matter 
power spectrum (in the Kaiser formula at a fixed scale, 
the non-Gaussian effects are perfectly degenerate with 
the galaxy bias). Furthermore, the redshift-space power 
spectrum of matter itself receives corrections of order 
/nl from non-Gaussianity, which are not present in real 
space. 

The non-Gaussian corrections to the Kaiser formula 
lead to the interesting question of whether there are de- 
generacies of the effects of non-Gaussianity with other pa- 
rameters measured from redshift-space distortions, such 
as dark ene rgy parameters or consistency tests of General 
Relativity [2a, [23| . The severity of the degeneracies will 
depend on the level of allowed non-Gaussianity. While 
the local model is likely to be constrained tightly in the 
near future, other models such as the equilateral model 
are much less constrained, but can still be lead to notice- 
able effects in redshift space (Fig. S]). 

Finally, we expect that non- Gaussian effects on halo 
velocities will have an even larger impact on the bispec- 
trum, of halos in redshift space. Again, these questions 
deserve more study via higher order perturbation theory 
as well as N-body simulations. 
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Appendix A: Matter Bispectra from Primordial Non-Gaussianity 

Primordial non-Gaussianity is most easily characterized by the bispectrum of the primordial curvature perturba- 
tions. Different shapes, i.e. configuration dependences, have been proposed in the literature. Since we work in the 
sub-horizon regime at late times, we phrase non-Gaussianity in terms of the potential $ at early times, which can be 
directly related to the primordial perturbations. The bispectrum B^ at a fixed redshift is defined by 

($(ki)$(k2)$(k3))=(2^)35,,(ki+k2+k3)S$(fci,/c2,A:3). (Al) 

Throughout, we take i?$ to be defined at last scattering ("CMB convention", z* w 1100). The bispectrum corre- 
sponding to local non-Gaussianity of the type Eq. ([3]) is 

B'^'{h,k2, h) =. 2/NL(P*(fci)P*(fc2) + 2 perm.), (A2) 

where P$ is the power spectrum of $ at z* . The bispectra corresponding to equilateral and (en) folded types are given 
by 

Bl""' = 6/nl \-P1P2 - 2 perm. - 2{PiP2Pzf'^ + P^^^P^^^Ps + 5 perm.l (A3) 



2/3 _ pl/3p2/3. 



BT = 6/nl [P1P2 + 2 perm. + 3(PiP2P3)'^' - P^ P2 P^ - 5 Perm.J . (A4) 

(A5) 

In order to convert this to the matter bispectrum smoothed with the window function WR^ we use the Poisson equation 
(again, valid on sub- horizon scales): 

3 (1 + z)H^\lm 

where T{k) is the matter transfer function, and ^^(z) oc {\-\- z)D[z) is the potential growth factor normalized to unity 
at last scattering. Then, power spectra and bispectra are related by 

PR{k,z) = Ml{k,z)P^{k), (A7) 

BR{ki,k2,k3;z) = MR{ki,z)MR{k2,z)MR{k3,z)B^{ki,k2,k3). (A8) 

Appendix B: Derivation of the velocity difference distribution 

In this appendix we present some details on the derivation of the relative velocity distribution for peaks, starting 
from Eq. (|29p . Expressing the Dirac and Heavyside functions in terms of their Fourier transforms, using 6r = —'Vur, 
and using Eq. p3p as well as partial integration, we can write the velocity probability distribution of peaks as: 

P(uri,Mr2;xi,X2) = / dai (ia2 rQ / — -^ exp(-za,.0'~) Z[J], (Bl) 

Jua-R J van \r=l'' J 

r(x) = 0ia'M^fl(|x-xi|) + 02a'Tyfl(|x-X2|) (B2) 

-H(/.3l^fl(|x - xi|)f' + 04M/fl(|x - X2|)f\ (B3) 

Here, we have defined 013 = u^i, 0^4 = Ur2 in order to allow for more compact notation, and used the Einstein summing 
convention. In the following we adopt the notation of [l3| . Using Eq. p4p and extracting the terms proportional to 
(j)^^ r = 1...4, leaves the (j) integrals as Fourier transforms of Gauss exponentials. We obtain: 



P(Wrl,Mr2;Xi,X2) = ,„ .„ „ / dui / da2 (B4) 



4 ^2 



1 4 ^2 oc z 1 VI ■* 



... ^ .^^^ SXj 

%^j -^ n— 3 
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Here, we have denoted ordered sequences of integers ii, . . . ,i„ with [in], and the sum runs over aU those ordered 
sequences with each index running from 1 to 4. Furthermore, we have defined 

4:l....„ ^ |n/d'x,W^(^;)(|x,-r,J)l ei:L„(xi,...x„), (B6) 



rk^^, _ J d'^Wni^), 1 = 1,2 



^«w - ..x.r:..:' . o\ (B7) 



_ (TH, i= 1,2 
^«,« = ^ _ -■ ^ 3 4 (B8) 

This is a formidable array of definitions, which however allows us to keep the treatment general. In particular, the 
extension to higher-order velocity distributions, or mixed peak density and velocity distributions is straightforward. 
In evaluating the derivatives with respect to aj , it is useful to define the following function: 

= -=L=H„,_i(i./V2)cxp('-iz.2\ (B^Q) 

where H„i denotes the Hermite polynomials. We then define the following coefficients: 

a™ = /„H«z/"-iexp(-z/V2), rn > I (Bll) 

ao = ^/^CTic{v/V2)^ly^'^cxp{-:y^/2) (B12) 

bm = /m+l(Wrl/cr„), m>0 (B13) 

Cm = Im+l{Ur2/<yu), TO > 0. (B14) 

Note that the a„i differ by factors of (27r)^/^ from those defined in [1J|. The second approximate equality in Eqs. ()Blip - 
(jB12p is valid in the high-peak limit, v ':$> \, which we will assume throughout. Of course, we do not make this 
approximation for the velocities Wri,'Sr2. The first few velocity coefficients are 

60 = exp(u2^/2a2) (B15) 

61 - — exp(u2^/2a2) (B16) 



a,, 



ul 



u. 



u 



rl_ 1 I „„„,'„- 2 /r,_2 
2 



- 1 exp(w^i/2a^), (B17) 



and correspondingly for cq, Ci, C2- 

For the peak correlation functions it is possible, by careful definition of different sets of integers, to expand the 
exponential in Eq. (|B4p , and reorder the terms into a sum over powers of a^ [l3| . In our case, this expression would 
become quite cumbersome however due to the different types of a parameters for peak density and velocity. Such an 
expression does not appear to be particularly useful. Instead, we make the assumption that the smoothed correlation 
functions Wj^ are much less than unity, so that only the linear term of the expansion of the exponential needs to be 
kept. This will be appropriate in particular in our case since we are interested in large scales where the correlations are 
small. It is then straightforward to perform the a integrals. The end result is that each derivative d"^ /dd^ (including 
n = 0) is replaced by a„, if i = 1, 2, 6„ if i = 3, and c„ if i = 4. 

Analogous to (but much simpler than) the above derivation, one can obtain an expression for the peak density 
Eq. (j27p in terms of the density ri-point functions: 

n^^{y) = (27r)-^/^ / daexp <( ^ ^-^<^(0, . . . , 0)^ 1 e'^^^ (g^g) 

(B19) 
(B20) 
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where for the second equality we have again assumed small correlations and truncated the series at the three-point 
function ^)j . In the third equality, we have further used the v ^ 1 limit. 

We can now assemble the velocity distribution 7^(uri,Mr2;xi,X2), collecting all non- vanishing two-point and three- 
point terms. First, we order terms in the following way: 



V{Url,Ur2;y^l,X2) = < A + B ^ ^^ 



Url ( U,.2 



D 






- 1 



Ur2 / Mrl 



- 1 



27r (7- 



■ exp 



The coefficients of the different velocity terms arc given by: 

■f 



A = i + hiei{T)+hi 



t5&& 



?r(o)+er(xi,xi,x2) 






^^ier (0) 



bL. 
cr,. 



2(Ju 



B = ^ei"M + :^[Ci'"(xi,xi,X2) + 2ei*«(xi,X2,X2)] 



c 

D 



5 1 5-?/? (,Xl,Xi,X2J, 

^^ u ^^ u 

-2^er"(xi,xi,x2), 



E = -|[^i""(0)-fe|r«(^l'^2,X2) 



(B21) 
(B22) 
(B23) 

(B24) 
(B25) 
(B26) 
(B27) 
(B28) 



Here, we have defined the first order Lagrangian bias b^ = v I ctr- The last term for A comes from the normalization 
by the number of peaks, Eq. (|B20[) . We can now evaluate Eq. (pO)) to obtain the peak correlation function (in real 
space), 

1 + epk(r) = 6ie^^(r) + bie,^\r), (B29) 

recovering the well-known expression from e.g. [g. Correspondingly, the real-space peak power spectrum reads 

Ppk(fc) - hlPnik) + bl/'\k). (B30) 

In order to obtain the matter- weighted velocities, we have to divide by (1 -I- ^pk) [Eq. (|3T|) ]. This yields 

A' = A- bli^{r) - 6i^i*'*'(xi, XI, X2) = 1. (B31) 

In this way, Vsu is properly normalized to second order in the correlations. Performing the integral in Eq. pip term 
by term, we obtain: 



P5«(<5m;xi,X2) = \l + B— +{E-C)-[--^-2 

Su f Su^ 






1 fSw 
4 \ at, 



4'T„ V <^l 
1 f Su^ 



(B32) 



Clearly, this expansion is equivalent to the Edgeworth expansion [Eq. (|^T|) ]. and the coefficients B through E are 
easily related to the moments of Su. 



Using Eq. (|16l) again to express the various three-point correlations, we have 

{6u) = 2Ba„ = 26iei"(r)+6iK^*"(r)-2C^"*(r)] 

(Su^) = 2il + E- C)al = 2[al + Cd'ir) + bLCn'iO) + bLCn^r) - 26iei™(r) 

{6u') = 12 I? a^ 

((5m*) - 3{5uy = 0. 



6a™M 
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(B33) 

(B34) 
(B35) 
(B36) 



These moments agree precisely with those derived in § IIII Al 
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